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Abstract. The hypothesis that dark matter and dark energy are unified through the 
Chaplygin gas, an exotic fluid obeying p = — A/p, is reexamined. Using generalizations 
of the spherical model which incorporate effects of the acoustic horizon we show that 
an initially perturbative Chaplygin gas evolves into a mixed system containing cold 
dark matter like gravitational condensate. 
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1. Introduction 

It has long been appreciated tliat there exists a dark-matter problem in that the 
nonrelativistic matter fraction Qm of critical density, inferred from galactic rotation 
curves, cluster stability, and peculiar velocities, far exceeds the ordinary baryonic 
matter fraction Qb — 0.04 provided by nucleosynthesis, but falls short of the standard 
inflationary prediction of unity pp. Observations of high-redshift supernova and 
microwave background anisotropies [SI E] have now settled the case for a flat universe 
with the dominant component dubbed 'dark energy' possessing appreciable negative 
pressure, such that the ratio w = p/ p is less than —0.8. The most conservative reading 
of the WMAP results |lj is that ^cdm — 0.23 is provided by traditional cold dark 
matter (CDM) candidates, while the balance resides in a cosmological constant A with 
the corresponding density fraction — 0.73. The inherent fine-tuning and coincidence 
problems of the ACDM model are somewhat ameliorated in quintessence models 
which replace A by an evolving scalar field. However, like its predecessor, a quintessence- 
CDM model assumes that dark matter and dark energy are distinct entities. For a recent 
review on the most popular dark- matter and dark-energy models, see jH]. 

A more liberal interpretation is that these data may tell us something profound 
about the space-time we inhabit, dark matter/energy being different manifestations of a 
common structure encoded in the dark sector's w [Zj. An interesting possibihty becomes 
evident upon recalling that the average density in a virialized halo is more than several 
hundred times higher than the critical density, suggesting a density dependent w = w{p) 
with w ^ for p ^ Per, ~ — 1 for p ~ Per and the identification of "dark matter" 
with objects which became nonlinear at high redshift. The general class of models in 
which a unification of dark matter and dark energy is achieved through a single entity 
is often referred to as quartessence |Hlin]- Among other scenarios of unification that 
have recently been suggested, interesting attempts are based on the so-called k- essence 
[TUl ITT] , a scalar field with noncanonical kinetic terms which was first introduced as a 
model for inflation [T2| . 

Perhaps the simplest scenario in which dark matter and dark energy are different 
manifestations of a single substance may be realized through the Chaplygin gas, an 
exotic fluid obeying 

which has been extensively studied for its mathematical properties The 
cosmological potential of equation was first noted by Kamenshchik et al Jlj who 
observed that integrating the energy conservation equation in a homogeneous model led 
to 



where a is the scale factor normalized to unity today and B an integration constant. 
Thus, the Chaplygin gas interpolates between matter, p ~ a~^, at high redshift and a 




(2) 
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cosmological constant like p ~ -\/]4 as a tends to infinity. 

Of particular interest is tliat equation (|T)) is obtained [121 [13 [HI from the 
Lagrangian 

c = -Va ^i-g^^-e,^ 9,, , (3) 

which may be seen by evaluating the stress-energy tensor T^^ and introducing = 
0,fj. I \J Q'^^ 6ifi for the four-velocity and p = \/~A/ ^1 — g'^'^ 6,^ 6,^ for the energy 
density. One recognizes £ as a Lagrangian of the Born-Infeld type familiar in the D- 
brane constructions of string/M theory j^. Geometrically, £ describes space-time as a 
3+1 brane in a 4-(-l bulk via the embedding coordinate = 9. The same Lagrangian 
appears as the leading term in Sundrum's [T7j effective field theory approach to large 
extra dimensions; for one such extra dimension, < 6* < £5, lowering the Planck scale 
to the intermediate scale of about 10^^ GeV, one estimates 

^1/8 

~ £5 ^ ~ meV jl6j . 

The Born-Infeld Lagrangian Q is a special case of the string-theory inspired 
tachyon Lagrangian ^H] in which the constant \fA is replaced by a potential V{9). 
Unfortunately, tachyon models encounter the problem that, generically, V{Q) is a 
maximum, so the tachyon field rolls quickly towards the minimum V^(oo) = thereby 
driving tf = 6'^ — 1 — > to leave only a dust component JH]- However, the situation 
may be quite different in inhomogeneous tachyon cosmology [201 • 

It is possible to define a 'generalized Chaplygin gas' [21] through 

V = -^/P", (4) 

where a > 0. In this case, the analogue of equation Q lacks any geometrical 
interpretation, but equation (0} can be obtained from a moving brane in Schwarzschild- 
anti-de-Sitter bulk [22] ■ In the limit a ^ 0, one approaches a model equivalent to ACDM 
in the background equations, although the perturbation equations yield a different power 
spectrum [23] • 

All models that unify dark matter and dark energy face the problem of nonvanishing 
sound speed and the well-known Jeans instability. A fluid with a nonzero sound speed 
has a characteristic scale below which the pressure effectively opposes gravity. Hence the 
perturbations of the scale smaller than the sonic horizon will be prevented from growing. 
In the Chaplygin-gas model, the sound speed is small at early times but is significant at 
the crossover region and hence the small scale perturbations are dumped. Indeed, it has 
been demonstrated that the purely perturbative Chaplygin-gas model [22[2S1 is in severe 
conflict with the cosmic microwave background (CMB) angular power spectrum [2^1 [2Z1 
and the mass power spectrum '28^ However, it has recently been pointed out that if 
nonadiabatic perturbations are allowed, the Chaplygin-gas model may be compatible 
with the 2df mass power spectrum p,9\ . 

A more general situation where the Chaplygin gas is mixed with CDM is considered 
in a number of papers [SHI IHH [H21 IHHl [211 [Ml [HH [HZ] • In this case, the Chaplygin gas 
simply plays the role of dark energy. In keeping with the quartessence philosophy it 
would be preferred if cold dark matter could be replaced by droplets of Chaplygin 
condensate. Homogeneous world models containing a mixture of cold dark matter 
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and Chaplygin gas have been successfully confronted with lensing statistics jHUJ EI] 
as well as with supernova and other tests [321 EH]- Bean and Dore and similarly 
Amendola et al [HSj have examined a mixture of CDM and the generalized Chaplygin 
gas against supernova, large-scale structure, and CMB constraints and demonstrated 
that a thorough likelihood analysis favors the limit a —>■ 0, i.e., the equivalent ACDM 
model. In both papers it has been concluded that the standard Chaplygin gas is ruled 
out as a candidate for dark energy. However, a recent analysis IHH] of the supernova 
data seems to indicate that the generalized Chaplygin gas with a > 1 is favored over 
the a — s> model. 

In spite of these somewhat discouraging results we share the opinion of Gorini et 
al inni that the Chaplygin gas still deserves further investigation and that a thorough 
investigation of the nonlinear regime of the growth of inhomogeneities is needed for a 
definite conclusion concerning the compatibility of the Chaplygin gas cosmologies with 
the observable large-scale structure of the universe. On the basis of an extended spherical 
model which includes the sonic horizon effects, in this paper we show that a fraction of 
the Chaplygin gas condensates and never reaches a stage where its properties change 
from dark-matter-like to dark-energy-like. Hence, at the crossover region a fraction of 
the gas behaves effectively as dark matter and the rest as dark energy. 

We organize the paper as follows. In section 2 we discuss the nonlinear evolution 
of perturbation in the Newtonian limit by generalizing the spherical model. Summary 
and conclusions are given in section 3. 

2. Condensate Formation 

To be able to claim that the Chaplygin gas (or any other candidate) actually achieves 
unification, one must be assured that initial perturbations can evolve into a deeply 
nonlinear regime to form a gravitational condensate of superparticles that can play the 
role of cold dark matter. In ITH] this was inferred on the basis of the Zel'dovich 
approximation [lU]. In comoving coordinates, the solution for the inhomogeneous 
Chaplygin-gas cosmology is 



Here 7 is the determinant of the induced spatial metric tensor •jij = giogjo/goo ^ 9ij and 
B can be taken as constant on the scales of interest. The generalization ^ of equation 
(0) allows us to implement the geometric version of the Zel'dovich approximation: the 
transformation from Euler to Lagrange (comoving) coordinates induces as 



where Di' is the deformation tensor, is the velocity potential, and the quantity h = b{t) 
describes the evolution of the perturbation. The Zel'dovich approximation offers a means 
of extrapolation into the nonperturbative regime via equation (0) and 




(5) 





(6) 




(7) 
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where the A, are the eigenvalues of ip/ . When one (or more) of the A's is (are) positive, 
a caustic forms on which 7 — > and — > 0, i.e., at the locations where structure forms 
the Chaplygin gas behaves as dark matter. Conversely, when all of the A's are negative, 
a void forms, p is driven to its limiting value \/A, and the Chaplygin gas behaves as 
dark energy, driving accelerated expansion. 

For the issue at hand, the usual Zel'dovich approximation has the shortcoming 
that the effects of finite sound speed are neglected. Since the structure formation occurs 
in the decelerating phase, we can address the question in the Newtonian limit p <^ p 
by generalizing the spherical model. In the case of vanishing shear and rotation, the 
continuity and Euler-Poisson equations become 

p + 3np = 0, (8) 

3W + 3H= + 4.Gp + |-gAp)=0. (9) 

where Cg is the adiabatic speed of sound given by 
2 dp A 

= TT- = (10) 



7i is the local Hubble parameter, and the time derivative is at the fixed Lagrangian 
coordinate r. It is convenient to introduce the comoving coordinates x = f/a. Then 
writing 

p = p{l + 6), (11) 

H = H + 6H, (12) 

subtracting the background, eliminating 60. and 57i, and changing the variables from t 
to a we find 

2../ 3 3,,^ 4{aST 1 + S d ( cl dS\ ^ 

a^6" + 2 - 2 + ^) - 3 ITT - (iT^ ^ j = °' ^''^ 

where ' denotes the derivative with respect to a. To solve this partial differential equation 
as it stands, one would need a rather involved numerical computation. Instead, we 
proceed first by discussing the linear approximation of ()13|) and then by solving the 
nonlinear problem imposing certain restrictions on the solution. 

First, consider the linear approximation which has been discussed in detail by Fabris 
et al Keeping only the terms linear in 5, 5', and 5" and using the expansion 

5(a,f) =^(5pert(fc,a)e^'=^ (14) 

k 

one obtains an explicit solution for the perturbative density contrast which may be 
expressed as 

5^en{k,a) oc a^^^'^J^/u^dsk). (15) 
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Here Ju{z) is the Bessel function, k the comoving wave number, and ds the comoving 
sonic horizon given by 



with Cs = V A/p and the equivalent matter fraction 



n=^B/iA + B) = ^/B/p,,. (17) 

Thus, for dgk <^ 1, Spen ~ a, but for dgk ^ 1, 5port undergoes damped oscillations. 
Similar results have been obtained numerically in the relativistic case |42j . 

Next, consider the nonlinear evolution. To simplify equation ()13p. it seems 
reasonable to make the following restriction: we take 6{t, x) to be of fixed shape, wave 
number dependent size, and time dependent amplitude, then evaluate the comoving 
coordinate derivatives at the origin. That is to say, we retain the spirit of the spherical 
dust model where each comoving coordinate point is treated as an independent origin, 
but allow for the local coupling of the inhomogeneity curvature to sound waves. A 
convenient choice is the spherical lump 

S{a,x)=SR{a)fix/R), (18) 

where f{z) is an arbitrary function satisfying /'(O) = 1 and /"(O) < 0. For simplicity, 
we take f{z) = exp(— 2:^/2). Evaluating (fT^ at the origin, we obtain an ordinary 
differential equation for 6ji 

a'5n + ^a6^--6Ril + 5R)--'^^ + - (-) j^^^ = 0, (19) 

where 



2^2"^ ' 3 1 + 5k 4 \anJ (l + 6r^ 



al' = {dM/Rf"a-^ = 8.06 (l - n^^fn-^' (R-^h-^ M^cf . (20) 
The recommendations of equation ()19|) are that it 

(i) reproduces equation (fT3j) at linear order and 

(ii) extends the spherical dust model by incorporating the Jeans length through the 
last term. 

To gain some perspective, with WMAP fixing Oeq = 3 ■ 10^^ and reionization 
at a redshift of 20, initial perturbations must satisfy 5(aoq) > 6 ■ 10~^ on mass 
scales M< W^Mq even for ACDM. Taking Vt = 0.27 and h = 0.71, we find a;^^ = 
15.4(_R"^ Mpc)^/^, which coincides with 1 + Zreion = 21 at i? = 0.338 Mpc corresponding 
to a protogalaxy mass M = 2.3 -IO^^Mq. 

In figure Q we show the evolution of two initial perturbations from the radiation- 
matter equality for or = 1/21. In contrast to linear theory, where for any R the 
acoustic horizon will eventually stop 6r from growing irrespective of the initial value 
of the perturbation, here, for initial SR^adec) above a certain threshold, Sr^o) — oo at 
finite a just as in the dust model. Conversely, at sufficiently small ^^('^dec), the acoustic 
horizon can stop 6R{a) from growing even in a mildly nonlinear regime. Figure El shows 
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10"^ 10"^ 10"^ a 

Figure 1. Evolution of SR{a) in the spherical model, equation (|19|l . from acq = 
3 X lO-'^ for an = 0.0476, 5ij(aeq) =0.004 (solid) and ^^(aeq) =0.005 (dashed). 

how the threshold ^^(adec) divides the two regimes depending on the comoving scale R. 
Qualitatively similar conclusions have been reached in a different way by Avelino et al 

m- 

One question is how accurate is the approximation (|TH|l which assumes that the 
comoving size is time independent? To check this, we have relaxed the assumption 
(fTKjl by allowing the second derivative of / at the origin to depend on a, i.e., the initial 
Gaussian perturbation to have a variable width. Evolution of the width is then described 
by a second-order differential equation coupled to equation |T9|l . By solving this system 
of coupled equations we find a slight change of the behavior of 6r in the linear regime 
but no significant change of the linear/nonlinear boundary depicted in figure |21 

The crucial question is what fraction of Chaplygin gas goes into condensate. Indeed, 
in jUj we have noted that if this were 92% as given by the geometric Zel'dovich 
approximation, the CMB and the mass power spectrum would be reproduced. To answer 
this quantitatively, we adopt the Press-Schechter approach Assuming Sji^adec) is 
given by a Gaussian random field with dispersion a"(i?), and including the notorious 
factor of 2 to account for the cloud in cloud problem, the condensate fraction at a scale 
R is given by 

F(R) = 2r expf-4l^Uerfc(^V (21) 
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10"^ 10"' 10' 10^ 



R[Mpc J 



Figure 2. Initial value 5fl(adcc) versus R for ^l = 0.27 and h = 0.71. The threshold 
^c(adoc) is shown by the line separating the two regimes. The dashed line gives cr(i?) 
calculated using the concordance model. 



where 6c{R) is the threshold (5c(adec) depicted in figure El In figure El we also exhibit the 
dispersion 

poo rjh 

a\R) = / — eM-k''R')^\k,aA,,), (22) 
JO k 

calculated using / above as the window function and the concordance model |4J variance 
A^(M)^co„st(A)V^w(^)"". (23) 

In figure 121 we present F{R) calculated using (|7T |) -(PH |) with const=7.11x 10~^, the 
spectral index n=1.02, and the parametrization of Bardeen et al KHl for the transfer 
function T{k) with r2B=0.04. The parameters are fixed by fitting (PHj) to the 2dFGRS 
power spectrum data jlTj. Our result, in qualitative agreement with that of Avelino et 
al jlS], demonstrates that the collapse fraction is less than 1% , far too small to affect 
the conclusions of |2Zl EHl • 



3. Summary and Conclusions 

Thus we have found that nonlinear condensate, while present, is insufficient to save the 
simple Chaplygin gas. In effect, the model is a victim of the radiation dominated phase 
which turns the Harrison-Zel'dovich spectrum 5k ~ k^^"^ to 5k ~ k~'^f'^ at i?eq — 26Mpc. 
In the pure Chaplygin-gas universe there would inevitably be insufficient small scale 
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Figure 3. Fraction of Cliaplygin gas in collapsed objects using Sc{R) and (j{R) from 
figure 121 



power to drive condensation. This is a more definitive, if pessimistic, statement tlian 
[43] because we have actually followed the perturbations into the fully nonlinear regime. 

Pragmatically, one can follow [Zj by introducing entropy perturbations to make the 
effective speed of sound vanish, even in the nonperturbative regime [23 • Caution should, 
however, be taken that imposing 5p = as an initial condition does not guarantee its 
propagation except in the long-wavelength limit |1H] where there is no problem in any 
case. More important is that the negative results were drawn without exploiting the 
braneworld connection of equation (jHl): in braneworld models (for a recent review, see 
|49j ) the Einstein equations are modified, e.g., by dark radiation and hence so are 
equations ()8|9|) . From a different perspective, similar changes are brought about by the 
radion mode which gives a scalar-tensor gravity. 

An interesting idea how to circumvent the difficulties with structure formation in 
the Chaplygin gas model has recently been put forward by Bento et al Following 
PU] , they have identified p = —p\ and pdm = P + P, so that the energy- momentum 
tensor appears as a combination of dark matter pdm and the cosmological "constant" 
A. Still, T^^ = supplies the continuity and Euler equations while p + 3p = pdm — 2pA 
sources the gravitational potential. In particular, {dp^jdri) / remains in the Euler 
equation. Since |p| < p ~ pd 

ni in the Newtonian context, and Cdj-Q — -dpA/dpdm = 
c^/ (1 + c^) ~ c^, there is no modification to equations (jHl) and Q, or the conclusions 
on structure formation. The reason Bento et al [HII conclude this bookkeeping allows 
^dm to grow unimpeded is that they have taken over the equations of Arcuri and Waga 
[52] . wherein A is assumed time dependent but homogeneous - c/ equation (14) therein. 
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Here, instead, A homogeneous implies pdm homogeneous, i.e. S^m = 0. 

Our prognosis is then that the crisis in quartessence cosmology is akin to that which 
existed in electroweak gauge theory before the Higgs mechanism. The Chaplygin-gas 
model is not so much wrong as incomplete. 
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